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O ' Abstract 

Q_ij We point out that thermal relic abundance of the dark matter is strongly 

Qh! altered by a non-perturbative effect called the Sommerfeld enhancement, when 

• constituent particles of the dark matter are non-singlet under the SU(2)i gauge 

interaction and much heavier than the weak gauge bosons. Typical candidates 

- T— I | 

^ ■ for such dark matter particles are the heavy wino- and higgsino-like neutrali- 

se 1 

nos. We investigate the non-perturbative effect on the relic abundance of dark 
matter for the wino-like neutralino as an example. We show that its ther- 
mal abundance is reduced by 50% compared to the perturbative result. The 
wino-like neutralino mass consistent with the observed dark matter abundance 
turns out to be 2.7 TeV < m < 3.0 TeV. 



1 Introduction 



Precise measurements of cosmological parameters have achieved amazing progress 
in recent years. In particular, the observation of cosmic microwave background 
anisotropics by the Wilkinson Microwave Anisotropy Probe (WMAP) confirmed 
that non-baryonic dark matter amounts to 20% of the energy of our universe. The 
existence of dark matter forces us to consider physics beyond the standard model 
(SM) for its constituent, because the SM has no candidate for the dark matter. 

Many models beyond the SM have been proposed for providing dark matter 
candidates. Models involving weakly interacting massive particles (WIMPs), such as 
supersymmetric models and universal extra dimension (UED) models [21 El; have an 
advantage over other models, since the WIMPs can explain the observed abundance 
naturally in the thermal relic scenario jSl S El • Comparison of the predicted thermal 
abundance with observations is a powerful tool to constrain those models. It is now 
important for discussion of new physics signatures at collider experiments such as 
the LHC. 

In the thermal relic scenario, the annihilation cross sections of dark matter play 
a central role in evaluation of the abundance. Perturbative calculation of the cross 
sections is sufficient for it in usual cases. However, a non-perturbative effect on 
the cross sections has to be included, when dark matter particles are non-singlet 
under the SU(2)^ interaction and much heavier than the weak gauge bosons. The 
weak interaction is not a short-distance force, but rather a long-distance one for 
non-relativistic two-bodies states of such heavy particles. The wave functions for 
the two-bodes states are modified from the plane waves by the interaction, and the 
annihilation cross sections are affected. The dark matter particles are non-relativistic 
at the freeze-out temperature. The attractive channels due to the weak force enhance 
the annihilation cross sections, compared with the perturbative ones. This effect is 
called the Sommerfeld enhancement, which was found in inelastic reactions between 
non-relativistic charged particles, historically. 

In this letter, we point out that the non-perturbative effect strongly alters the 
relic abundance of dark matter compared to that in perturbative calculation. We 
investigate the non-perturbative effect quantitatively in the wino-like neutralino dark 
matter as an example, and show that the abundance is reduced by 50% due to the 
effect. We also discuss non-perturbative effects on other dark matter candidates at 
the end of this letter. 
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2 Sommerfeld enhancement in wino-like neutralino 
dark matter annihilation 



Winos are the SU(2)^ gauginos in the supersymmetric standard model, and they are 
an SU(2)i triplet. They are mixed with higgsinos, superpartners of the Higgs bosons, 
after the electroweak symmetry breaking. The neutral components in the mass 
eigenstates are called neutralino. The mixing with higgsinos is small for heavy winos, 
since it is suppressed by the electroweak scale. Thus, the wino-like neutralino is 
highly degenerate with its charged SU(2)l partner (x~) in mass when their masses 
are heavy enough. For the wino-like neutralino mass (m) of the order of 1 TeV, 
the mass difference between the neutralino and its SU(2)^ partner is dominated by 
a radiative correction [B]. It is calculated as 5m ~ 0.17 GeV, which is used for 
numerical calculations in this letter. 

There are four processes related to the calculation of the wino-like neutralino relic 
abundance, and those are X~X + 1 X°X~> an d X~X~ annihilation. We assume 

the CP conservation, and cross sections of x°X + an d X + X + annihilation are the same 
as those of x°X~ an d X~X~ > respectively. Furthermore, each process is decomposed 
into two ones with 5 = and 5 = 1, where S is the total spin of the two-bodies 
system at the initial state. In x°X an d XT XT annihilation, the S = 1 processes are 
forbidden at the s-wave annihilation. After all, processes we have to consider are 
X°X°i X~X + 1 X°X~ > XTX~ annihilation with S = 0, and those of X~X + 1 X°XT with 
S = 1. 

While gamma rays from the wino-like neutralino annihilation in the Galactic 
center (HI Ej and anti-particles fluxes from that in Galactic Halo jH] are evaluated 
including non-perturbative effects, the relic abundance of the neutralino have been 
calculated only within a perturbative method. However, as we will see later, the non- 
perturbative effect can significantly alter cross sections relevant to the abundance. 

Two-bodies effective Lagrangian, describing relative motion between two particles 
in the two-bodies system, is useful to evaluate the non-perturbative effect in the 
annihilation. The Lagrangian has a following form, 



s z \ x , r ) 



V 2 V 2 
id x o + -^ + -^- V(r) + 2iT5(r) 
Am m 



®s z (x, r) 



where x is the center of mass coordinate for the two-bodies system, while the relative 
one is represented by r. The index Sz is for the z-th component of the total spin S. 
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The potential term V(r) describes forces acting between two particles in the system. 
The absorptive part F stands for the annihilation of the two-bodies system. For 
derivation of the effective Lagrangian and following evaluation of the annihilation 
cross section using it, see Ref. j^j. 

In the x~X~ annihilation process with 5 = 0, the potential V(r) and absorptive 
part T in Eq. (JTJ) turn out to be 

V{r) = * + a 2 cl e - , r = ^|, (2) 

r r 2m z 

where a is the fine structure constant, a 2 = g 2 /^ is for the SU(2) gauge coupling 
constant, irtz is the Z boson mass, and cw = cosOw is for the weak mixing angle 
8w The first term in the potential comes from the exchange of photons, while the 
second one is from the exchange of Z bosons. In the calculation of T, we consider 
only the final states of SM particles, and neglect their masses, since they are light 
enough compared to the wino-like neutralino we are discussing. 

In both cases of the x°X~ annihilation processes with 5 = and 1, V(r) is 
induced from the exchange of W bosons, and both cases have the same form. On 
the other hand, the absorptive part T is different each other. These are given by 

-it i \ e~ m ^ r lna 2 2 25™* 

V(r) = -a 2 ^-, r>= 0) = --^, 1^) = -— , (3) 

where my/ is the W boson mass. The potential and absorptive terms in the xTX + 
annihilation with S = 1 are 

V r = a2C i r = ——1 . (4) 

r r 24 m z 

The x°X two-bodies system is mixed with x~X + state with S — 0, in which 
mixing occurs through a W boson exchange. Thus, the potential and absorptive 
terms are written by 2 x 2 matrices as 



V(r) 



/ a e ~ m zr e -m w r \ 

28m a 2 Cm — y/2a 2 

e ~m w r 

y -V2a 2 —^— J 




(5) 



Off-diagonal elements describe the transition between x~X + an d X°X° systems. 

As seen in these potentials, all processes have attractive channels except for 
that of x~X~ ■ The overlap between the wave functions of the incident particles are 
increased compared to the case without including the potentials, and it leads to 
enhancement of the annihilation cross sections. 



3 



Once the two-bodies effective Lagrangian is obtained, annihilation cross sections 
including the non-perturbative effect can be calculated through the formula, 

av = cT\A\ 2 , A = Jd 3 r e"^ + G(f, 0) , (6) 

where c = 2 for an annihilation of identical particles, otherwise c = 1. The relative 
velocity between incoming particles is denoted by v. The Green function G satisfies 
the equation of motion of the effective Lagrangian, 



mv 2 V 2 



+ V(r) + 2iV5 s (r) 



G{r,r*) = 8 3 {r-r*) . (7) 



4 m 

The boundary condition for the Green function is determined from following condi- 
tions. First, the Green function is analytic at any f and f except the point r — r 1 . 
Second, only out-going waves survive at large \r — r*|. 

In x°x° and x~X + annihilation processes with S = 0, the Green function has a 
2x2 matrix form. For this case, it has been found in Ref. IZI that above formula 
is simply extended as 

0%V = A ij T jj' A ij> ■ ( 8 ) 

3,5' 

The factor \A\ 2 in Eq. (jSJ is called the Sommerfeld enhancement factor, and it 
includes the non-perturbative effect due to the long-distance forces in V(r). Note 
that if we neglect the non-perturbative effect, the enhancement factor becomes one 
and the annihilation cross section in Eq. (jUJ) coincides with the result in a usual 
perturbative method. 

3 Effect of Sommerfeld enhancement on dark mat- 
ter abundance 

Now we evaluate the thermal relic abundance of the wino-like neutralino dark matter, 
including the non-perturbative effect. In the evaluation we have to include coannihi- 
lation processes in addition to the wino-like neutralino pair annihilation. We use the 
method developed in Ref. [3 El for the calculation of the relic abundance including 
coannihilation effects. Under reasonable assumptions, the relic density obeys the 
following Boltzmann equation, 



dY (<y c s'v) f x dg* 



dx Hx \ 3g* s dx 

4 



) s (Y 2 - Y 2 q ) . 



The yield of the dark matter, Y, is defined as Y = n/s, where n is the sum of the 
number densities of x°, X ~i an d X + ■ The variable, x = m/T, is the scaled inverse 
temperature of the universe. The equilibrium abundance, Y eq , is given by 

y eq 4n 1/2 ^v, (io) 

where g e $ is the number of the effective degrees of freedom defined as 

g cS (x) = 2 + 4(1 + 5m/m) 3/2 e- x5m/m . (11) 

The entropy density s and the Hubble parameter H are given by 

2tt 2 m 3 fg*y/ 2 7t m 2 



45 x 3 ' V10/ 3M P1 x 2 

where Mpi = 2.4 x 10 18 GeV is the reduced Planck mass. The relativistic degrees 
of freedom of the thermal bath, g* and g* s , should be treated as a function of the 
temperature for deriving the correct dark matter abundance in our calculation. 

The most important quantity to determine the abundance is the thermally aver- 
aged effective annihilation cross section (a e gv) in Eq. (j^J), defined as 

(<Tesv) = ^(a^)^ T (l + A J ) 3/2 (l + A,) 3/2 exp[-x(A, + A J )] , 



/ Tfi \ 3/*^ f f TTIV 2 \ 

( a H v ) = {^f) J 47lv dv (° r ^) ex P ( — ) > ( 13 ) 

where i, j = x°, X~ an d X + ? = = ^ m / m > A^o = 0, and is the annihila- 
tion cross section between i and j. In Fig. Q the enhancement ratio of the averaged 
cross section, (cr e sv), to that in the perturbative calculation is shown as a function 
of m with fixed m/T = 20, 200, 2000 (left figure) and as a function of T with fixed 
m = 2.8 TeV (right figure). For comparison, the cross section in a perturbative cal- 
culation is also shown as a dotted line in the right figure. Note that the perturbative 
cross section is constant in time. The little drop at x ~ 10 5 is due to the decoupling 
of x^- m the calculation of the cross section, we used the running gauge coupling 
constant at the m and mz in the absorptive terms and the potentials, respectively. 

In these figures, large enhancement of the cross section is found due to the non- 
perturbative effect when m is larger than ~ 1 TeV. A significant enhancement is 
shown at m ~ 2.4 GeV. This originates in the bound state composed of x°X° an d 
XTX + pairs jSlEl- The enhancement by the non-perturbative effect is more signif- 
icant for lower temperature. Since x° an d X~ are more non-relativistic for lower 
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Figure 1: Mass dependence of averaged cross section, (<r e fr), normalized by the perturba- 
tive one, (cr e fi)Ttee, for m/T = 20, 200, 2000 (left figure), and temperature dependence 
with m = 2.8 TeV (right figure). The perturbative result is also shown as a dotted line for 
comparison. Here, mass difference between x° an d is & m — 0-17 GeV. 



temperature, the long-range force acting between these particles strongly modifies 
their wave functions and alters the cross section significantly. 

Since the averaged cross section depends on temperature in a non-trivial way 
as shown in Fig. Ql we should integrate the Boltzmann equation numerically. After 
calculating the present value of the yield, Y Q , by the integration, we obtain the 
dark matter mass fraction in the current universe through the relation QoAih 2 = 
msoYoh 2 /p c , where p c is the critical density, p c = 1.05 x 10~ 5 /i 2 GeVcm -3 (h = 
0. 73j^o3) [HI; an d sq is the entropy density of the present universe. 

The result is shown in Fig. |21 In the left figure, the ratio of the yield with 
the enhancement to one without the enhancement (perturbative result) is shown as 
a function of temperature. The mass of the wino-like neutralino is assumed m = 
2.8 TeV. The enhancement of the cross section at the departure from the equilibrium 
decreases the abundance by 20 — 30%, and it leads to quick deviation of the yield 
from the perturbative one after decoupling. Since the annihilation cross section is 
increased for lower temperature by the Sommerfeld enhancement, the sudden freeze- 
out phenomenon on the yield does not occur and the yield continues to be reduced by 
the non-perturbative effect even for x > 100 compared to the perturbative one. The 
resultant dark matter abundance is reduced by 50% compared to the perturbative 
result. 

In the right figure, the relic abundance of the dark matter in the present universe 
is shown as a function of m in terms of Qh 2 . The allowed regions by the WMAP at 1 
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Figure 2: Ratio of yield with the non-perturbative effect to that in the perturbative cal- 
culation (left figure). Wino-like neutralino mass is fixed 2.8 TeV. Thermal relic abundance 
of the dark matter in the current universe as a function of wino-like neutralino mass (right 
figure). Allowed regions by the WMAP at 1(2) a levels are also shown as the dark (light) 
shaded area. 

and 2 a are also shown as shaded areas in this figure. We found that the mass in the 
wino-like neutralino dark matter consistent with the observation is shifted by 600 
GeV due to the non-perturbative effect and the wino-like neutralino mass consistent 
with WMAP results turns out to be 2.7 TeV < m < 3.0 TeV. 

4 Summary and discussion 

In this letter, we have pointed out the thermal relic abundance of dark matter, which 
is SU(2) l non-singlet and has a much larger mass than that of the weak gauge bosons, 
can be strongly reduced by the non-perturbative effect. We have investigated the 
non-perturbative effect on the relic abundance of wino-like neutralino as an example. 
Compared with the perturbative result, this effect reduces the abundance by about 
50% and increases the mass of the wino-like neutralino dark matter consistent with 
the observation by about 600 GeV. As a result, the thermal relic abundance of 
the wino-like neutralino dark matter is consistent with observed abundances when 
2.7 TeV < m < 3.0 TeV. 

The non-perturbative effect can change relic abundances of other dark matter 
candidates with 811(2)^ charge and heavy mass, such as higgsino-like neutralino. The 
non-perturbative effect on the thermal relic abundance of higgsino-like neutralino is 
expected to be roughly 10%, since winos are triplet under the S\J(2) L gauge group, 
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while higgsinos are SU(2)l doublet. Therefore, non-perturbative effect may change 
the abundance by 0(10%) for other SU(2)i doublet candidates for the dark matter. 
A detailed analysis of this subject is studied elsewhere. 

The Sommerfeld enhancement occurs reasonably for particles with electric charge. 
In fact the non-perturbative effect through photon exchanges for charged particle an- 
nihilation can change the abundance by about 10%. Therefore, one may think that 
the relic abundance of dark matter will be changed in the stau coannihilation region 
or in the case that gravitino is the lightest supersymmetric particle and it is pro- 
duced through decay of stau. However, the non-perturbative effect for stau does not 
change the dark matter abundance. In stau annihilation, a(f + f~) ~ a(f + f + ) and 
the Sommerfeld enhancement is positive (negative) for the former (latter) process. 
Hence, the enhancement is almost canceled and the change of the abundance is 1% 
at most. 

Kaluza- Klein (KK) right-handed leptons (E^) in UED models are highly degen- 
erate with the lightest KK particle in mass and a(E^E^) > a(E^E^) j31 IT2]. 
Hence, the enhancement could be expected to change the abundance of the KK dark 
matter. The change of the abundance is within 4% since a(E^E^) contributes to 
the effective annihilation cross section by 40% at most. 

Finally, we comment on the non-perturbative effect of colored particles, such as 
gluino. The enhancement for colored particles are very effective However, it 

may be very complicated due to the existence of the QCD phase transition, which is 
discussed in Ref. jHj. Furthermore, colored particles can not be candidates for the 
dark matter and are not expected to be degenerate with a dark matter particle in 
mass due to large radiative corrections by the strong interaction. Hence, this subject 
is beyond the scope of this letter. 
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